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O ■ Abstract. 

I The structure of the crust of a neutron star is completely determined by the experimentally 

■ measured nuclear masses up to a density of the order of lO'' g.cm^^. At higher densities, the 
^ ' composition of the crust still remains uncertain, mainly due to the presence of "free" superfluid 
O . neutrons which affect the properties of the nuclear "clusters". After briefly reviewing calculations 

■ of the equilibrium structure of the crust, we point out that the current approach based on the Wigner- 
Seitz approximation does not properly describe the unbound neutrons. We have recently abandoned 
this approximation by applying the band theory of solids. We have shown that the dynamical 
properties of the free neutrons are strongly affected by the clusters by performing 3D calculations 
with Bloch boundary conditions. 
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INTRODUCTION 



At the end point of the stellar evolution, neutron stars are the compact remnants of core 
collapse supernova explosions. Bom with temperatures as high as 10^^ — 10^^ K, the star 
rapidly cools down by emitting neutrinos and photons. A few hours after its birth, the 
^ , temperature of the star falls below about 10^ K and the external layers crystallize into 

a solid crust (for a general review of neutron star crust, see Within hundreds of 

years, the interior of the star becomes isothermal with temperatures typically less than 
10^ K(~0.1 keV). 

rS I From the nuclear physics point of view, a neutron star is a huge nucleus containing 

c3 ' about A ~ 10 baryons and a proton fraction of the order of 10%. A rough estimate of 

the radius and the mass of the star from the liquid drop model yields R = roA^^^ ~ 10 km 
and M = Airip ~ Mq respectively. The solid crust surrounding the star plays more or less 
the role of the neutron skin in heavy nuclei. Indeed correlations have been established 
between the neutron skin in lead ^ospj, ^j^g density at which the crust melts into 
a uniform liquid [3]. This gross picture however should not be taken too far since the 
conditions prevailing inside neutron stars are very different from those inside isolated 
nuclei. 

The crust of a neutron star represents only a few percent of the mass of the star 
but plays a crucial role for its evolution (magnetic field, cooling, bursts, starquakes, 
spinning-down, glitches, free precession, non axial deformations giving rise to the emis- 
sion of gravitational waves). We shall briefly review calculations of the equilibrium 
structure of the crust in the first section. The next section will be devoted to the neu- 
tron superfluidity in the crust. In the last section, we will show how the description of 



the inner crust can be improved by applying the band theory of solids and we will discuss 
some recent results. 



STRUCTURE AND COMPOSITION OF THE CRUST 



In the following it will be assumed that the complete thermodynamical equilibrium at 
zero temperature with respect to all interactions has been reached and therefore the 
matter is in its lowest energy state (this excludes newly-born hot neutron stars and 
neutron stars accreting matter from a companion star which will not be discussed here). 
This assumption is usually known as the cold catalyzed matter hypothesis. The ground 
state of the crust is obtained by minimizing the total energy per nucleon under the 
assumption of beta equilibrium and electroneutroneutrality. The crust is further supposed 
to be formed of a perfect crystal with a single nuclear species at lattice sites (for a 
discussion of possible deviations from this idealized model, see for instance [4|] and 
references therein). 



At densities below ~ 10^ g.cm^, the ground state of matter is a mixture of electrons 
and iron ^^Fe (the atoms are fully ionized at densities above ~ 10^ g.cm^^). At higher 
densities, nuclei become increasingly neutron rich due to inverse beta decay. Following 
the classical paper of Baym, Pethick and Sutherland [|5|], the total energy density in a 
given layer can be written as 



where n^^ is the number density of nuclei, E{A,Z} is the energy of a nucleus with Z 
protons and A — Z neutrons, £e is the electron energy density and £l is the lattice energy 
density. At densities p ^ 10^ g.cm^^ the electrons can be described as a relativistic 
Fermi gas. Assuming point like nuclei (the lattice spacing being very large compared to 
the size of the nuclei), the lattice energy density can be expressed as 



where c is a coefficient which depends on the lattice structure. For cubic structures, this 
coefficient is respectively equal to 0.89593, 0.89588 and 0.88006 for body centered, 
face centered and simple cubic lattice which suggests that the crust crystallizes in a 
body centered cubic lattice. 

The main physical input is therefore the energy of a nucleus. The structure of the outer 
crust is completely determined by the experimental nuclear masses up to a density of the 
order p ~ 6 x 10^^ g.cm^^ [6]. At higher densities the nuclei are so neutron rich that the 
energy E{A,Z} must be extrapolated. The composition of the nuclei in these layers is 
thus model dependent. Nevertheless most models predict the existence of nuclei with the 
magic neutron numbers = 50, 82, thus revealing the crucial role played by shell effects. 



outer crust 



Etot = nNE{A, Z} + £e + El 



(1) 




(2) 



TABLE 1. Sequence of nuclei in the outer crust of non-accreting cold neu- 
tron stars calculated by Riister et al. |0] (a theoretical nuclear mass table was 
used for the lower part). The last line corresponds to the neutron drip point. 
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The nuclei present at the bottom of the outer crust may be experimentally studied in the 
near future by several facilities, such as for instance FAIR at GSlQ, IS AC at TRIUMPH 
I, SPIRAL 2 at GANIL 1 and by the RIA project 0. The structure of the outer crust is 
shown in tabled] for one particular representative recent model. 



inner crust 

At the bottom of the outer crust, the nuclei become so neutron rich that some neutrons 
are no longer bound. The nuclear lattice then coexists with a neutron gas. The transition 
occurs when the electron Fermi energy becomes comparable to the binding energy of the 
protons in the nuclei. At the drip, the neutron excess 5 = (1 — 2A/Z) can be estimated 
as [H 

5drip = x/l -£o/5o - 1 , (3) 

where Eq is the energy per nucleon of infinite symmetric nuclear matter and 5o the 
symmetry energy. To lowest order in 5, the density threshold for the onset of neutron 
drip is approximately given by 

Pdrip - (^O^drip)^ X 10'^ g.Cm-3 . (4) 
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FIGURE 1. Nucleon number densities (in fm^^) along the axis joining two adjacent Wigner-Seitz 
cells of the inner crust of neutron stars for a few baryon densities iii, (in cm^^) as calculated by 
Negele&Vautherin 



For £o = -16 MeV and 5o = 32 MeV, we find 5drip ^ 0.225 and pdrip ~ 4 x lO^i g.cm^^ 
which is in remarkable agreement with the value obtained from more realistic nuclear 
models. The simple estimates ([3]) and @) illustrate the importance of the symmetry 
energy on the physics of neutron stars (see also [7] and references therein). It has been 
very recently shown using a liquid drop model that the composition of the inner crust is 
very sensitive to the density dependence of the symmetry energy [8]. 

The inner crust of a neutron star is a unique environment which is not accessible in 
the laboratory due to the presence of the "free" neutron gas. In the following we shall 
thus refer to the "nuclei" in the inner crust as "clusters" in order to emphasize these 
peculiarities. The structure of the inner crust has been studied using various approaches, 
mainly liquid drop and semi-classical models (for a recent review, see for instance 

The most realistic calculations of the structure of the inner crust were pioneered 
by the work of Negele&Vautherin [9]. Until very recently flO^, it was the only fully 
self-consistent quantum calculation. Expanding the density matrix in relative and cm. 
coordinates, Negele&Vautherin derived a set of effective equations for the nucleons, 
which closely resemble those obtained in the Hartree-Fock approximation with Skyrme 
forces. They determined the structure of the inner crust by minimizing the total energy 
per nucleon in the Wigner-Seitz sphere, treating the electrons as a relativistic Fermi gas. 
As pertains the choice of boundary conditions, they imposed that wavefunctions with 
even parity (even £) and the radial derivatives of wavefunctions with odd parity (odd £) 
vanish on the sphere r = Rceii- This prescription yielded roughly uniform density outside 
the nuclear clusters. The remaining spurious fluctuations were removed at each iteration 
step by averaging the densities in the vicinity of the cell edge. 

The composition of the crust is shown on table [2l These results are qualitatively 
similar to those obtained with liquid drop and semiclassical models. The remarkable 
distinctive feature is the existence of strong proton shell effects with a predominance of 
nuclear clusters with Z = 40 and Z = 50. Neutron shell effects are also important (while 
not obvious from the table) as can be inferred from the density fluctuations inside the 
clusters on figure [B This figure also shows that shell effects disappear at high densities 
where the matter becomes nearly homogeneous. 



TABLE 2. Sequence of nuclei in the inner crust of non- 
accreting cold neutron stars calculated by Negele & Vautherin Jgt]. 
is the total number of neutrons in the W-S sphere. 
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"pasta" phases 

The equilibrium composition of the clusters is the result of the competition between 
Coulomb and surface energies. At low densities, the lattice energy ([21) is a small contri- 
bution to the total Coulomb energy and nuclei are spherical. However at the bottom of 
the crust, the size of the cluster is of the same order of magnitude as the lattice spacing 
and consequently the lattice energy represents a large reduction of the total Coulomb 
energy (this reduction is about 15 % at the neutron drip). This means that the nuclear 
clusters may be strongly deformed in the high density layers of the inner crust. Rea- 



soning by analogy with percolating networks, Qgasawara&Sato llllll suggested that a 
transion to an "infinite network of linked nuclei" might occur at the bottom of the crust. 
The possibility of non spherical nuclear clusters, referred as "pastas", was considered by 
Ravenhall et al. [12] and Hashimoto &Yamada [ [131 ] who found from compressible liquid 
drop models, that as the density increases, the nuclei merge into cylinders ("spaghetti") 
followed by slabs ("lasagna"), cylindrical tubes and bubbles ("swiss cheese"). The pasta 
phases cover a small range of densities near the crust-core interface. Nevertheless they 
may represents up to half of the mass of the crust. The existence of these phases may 
have important astrophysical consequences for the gravitational wave emission and for 
pulsar glitches by changing the elastic properties of the crust[14], for the cooling of 
neutron stars by allowing direct URCA processes [15] and enhancing the heat capacity 
[3 17] and for core-collapse supernovae [18]. These pasta phases have been studied by 



various nuclear models, from liquid drop calculations to quantum molecular dynamic 
simulations (for a recent review see for instance [19] and references therein). However 
a few models do not predict the existence of such phases [ii 0, 0, Q. The energy 



differences between the various shapes are very small, typically less than ~ keV.fm"^, 
and as a result the structure of the crust is very sensitive to small differences between 
nuclear models. 



NEUTRON SUPERFLUIDITY IN THE CRUST 



The possibility of superfluidity in neutron stars was suggested a long time ago by Migdal 
[|24|l . only two years after the formulation of the BCS theory of electron superconductiv- 
ity. Microscopic calculations performed in pure neutron matter indicate that at densities 
below saturation density, neutrons are bound in Cooper pairs which condense into a ^Sq 
superfluid phase. However the density range for superfluidity and the magnitude of the 
pairing gap still remain uncertain due to different approximations of medium polariza- 
tion and self-energy effects which tend to suppress the pairing as compared to mean field 
calculations with bare nucleon-nucleon forces (for a review see for instance [ 25Ll26|l '). 

The situation is even more uncertain in the crust owing to the small nuclear asymmetry 
and to the presence of inhomogeneities (for a recent review, see ^Z^A ). Superfluidity of 
the "free" neutrons induces superfluidity of the bound neutrons inside the clusters and 
vice versa. This proximity effect tends to smooth the spatial variations of the pairing field 
(see and references therein). The effects of neutron superfluidity on the equilibrium 
structure of the crust have been considered by Baldo et al. n2S\\ in the density functional 
theory generalized to account for nucleon pairing. The nuclear lattice was treated in 
the Wigner-Seitz approximation. They determined the ground state of the crust at the 
baryon density p ~ 1.9 x 10^-^ g.cm^^, for which the neutron pairing is expected to be 
the strongest. They found that the structure of the crust is significantly affected by the 
pairing. Indeed, the charge of the cluster Z ~ 52 and the radius of the Wigner-Seitz 
sphere 7?ceii — 32 fm are increased by about 5Z ~ 8 and 5i?ceii — 3.5 fm respectively 
compared to calculations without including the pairing. However they emphasized that 
the results are very sensitive to the choice of the energy functional. These conclusions 
have been very recently confirmed by calculations at other densities [10]. This urge the 
need for a better understanding of pairing correlations in inhomogeneous neutron star 
crust matter. 



FROM NUCLEAR TO SOLID STATE PHYSICS 
Wigner-Seitz approximation 

In the quantum calculations briefly reviewed in the previous sections, the nuclear 
lattice was treated in the Wigner-Seitz approximation. This approximation however does 
not properly takes into account the unbound neutrons which are artificially confined 
inside the Wigner-Seitz sphere. This approximation leads to spurious fluctuations of 
the neutron densities and of the neutron pairing field [2?]. As a result calculations of 
the equilibrium structure of the crust and of crustal superfluidity are contaminated by 
unphysical shell effects which are very sensitive to the choice of boundary conditions 
that are imposed on the sphere [fsoll . These shell effects of the order of h^/2mR^^^^ where 
i?ceii is the radius of the Wigner-Seitz sphere, may be very large in the deep layers of 
the crust where the nuclear clusters are very close to each other. A correct treatment of 
the "free" neutrons requires the application of the band theory originally proposed for 



describing electrons in solids and recently applied to neutron star crust yjj, 132|] 



FIGURE 2. Energy spectrum (in MeV) of unbound neutrons in the outermost layers of the inner crust 
of neutron stars for different symmetry directions in k-space. The dashed Hne is the energy spectrum of 
the ideal Fermi gas. See reference (32] for details. 



Band theory 

The nuclear lattice can be partitioned into identical Wigner-Seitz cells whose shape is 
determined by the geometry of the lattice. The boundary conditions on the cell are not 
arbitrary but are fixed by the Floquet-Bloch theorem. Each single particle quantum state 
is described by a discrete index a and by a wave vector k. The energy spectrum is thus 
formed of a series of sheets or "bands" in k-space as illustrated on figure |2l We have re- 
cently applied the band theory by performing 3D calculations with Bloch boundary con- 
ditions (see [[32I1 and references therein), in order to investigate the effects of the nuclear 
clusters on the dynamical properties of the neutron gas. For this purpose, first ignoring 
the effects of pairing, we have studied the topology of the Fermi surface by computing 
its area and comparing it to the area of the corresponding sphere (remembering that the 
Fermi volume only depends on the density). We found that at low densities, meaning that 
the Fermi wave length of the unbound neutrons is much larger than the lattice spacing, 
the Fermi surface is nearly spherical, which is analogous to alkali metals. However at 
higher densities the Fermi surface is strongly deformed owing to Bragg diffraction as it 
is also observed for transition metals. This entails a renormalisation of the neutron mass 



into an effective mass defined by = pn/vn and given by a Fermi surface integral [|33|] . 
This effective mass can take very large values, up to about ~ 15mn at the baryon den- 
sity rijy = 0.03 fm^-^ [31]. We have also shown that in nearly homogeneous neutron star 
crust matter for which the BCS approximation is valid, the effects of neutron pairing on 
the effective mass are small and vanish in the limit of a uniform system 0411 . However 
as pointed out recently by Magierski fs^, strong inhomogeneities may lead to neutron 
localisation, which might further enhance the effective neutron mass. 



CONCLUSION 



The evolution of a neutron star is intimately related to the properties of its solid crust. 
Up to a density of about 10^ g.cm^, the crust is formed of a body centered cubic crystal 
of iron ^^Fe, which are fully ionized above ~ lO'^ g.cm^-^. With increasing density, the 
nuclei immersed in a relativistic electron gas, become more and more neutron rich owing 
to electron capture. At densities around 10^^ g.cm^^ neutrons start to drip out of nuclei 
and may form Cooper pairs which condense into a superfluid phase. At the bottom of 
the crust at densities of the order of 10^^ g.cm^, some calculations predict that nuclei 
may adopt exotic non spherical shapes, referred as "pastas". 

Unlike the nuclei in the shallower layers, the nuclear "clusters" in the inner crust 
cannot be studied in the laboratory due to the presence of "free" superfluid neutrons. 
The structure of the crust has been studied using different approximations and nuclear 
models. The current state-of-the-art of self-consistent quantum calculations which were 
pioneered by Negele&Vautherin[9], is the Hartree-Fock-Bogoliubov approximation. It 
has been recently shown in this framework that the neutron superfluidity greatly affects 
the composition of the crust[28]. 

However much remains to be done. Indeed, in all these calculations the lattice is 
treated in the Wigner-Seitz approximation in which the dripped neutrons are artificially 
confined in spheres. A much more accurate description of the crust, which is essential 
in order to interpret observations of neutron stars, should rely on the band theory of 
solids. We have recently shown that the dynamical properties of the superfluid neutrons 
are strongly affected by the nuclear clusters by carrying out 3D calculations with Bloch 
boundary conditions 1321]. This work is a first step towards realistic calculations of the 
properties of neutron star crust. 
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